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^ ' Abstract 

Within the 5D standing wave braneworld model numerical solutions of the equations for 

matter fields with various spins are found. It is shown that corresponding action integrals 

'^ I are factorizable and convergent over the extra coordinate, i.e. 4D fields are localized on 

j~^ ■ the brane. We find that only left massless fermions are localized on the brane, while the 

^D , right fermions are localized in the bulk. We demonstrate also quantization of Kaluza-Klein 

excited modes in our model. 
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><! : 1 Introduction 



Braneworld models with large extra dimensions jHll] have attracted wide attention owing to 
their success in solving some open questions in modern physics J5HH]- To realize the braneworld 
idea one should find a mechanism of localization, however, there are some difficulties in this 
matter within standard brane approaches [9lll0). It is especially economical to consider models 
with a pure gravitational localization mechanism, since gravity has universal couplings with all 
matter fields. In this paper we consider one such scenario within the standing waves braneworld 
[TTiri2] and explicitly show existence of the universal gravitational trapping of all kinds of matter 
fields on the brane. Main differences of the model with the most brane scenarios are that our 
metric ansatz is time-dependent and contains increasing warp factor. 
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The standing waves braneworld is formed by the background solution [TT||T2]: 

ds^ = e^'^l^i [dt^ - e^dx^ - e'^dy^ - e-^^'dz^) - dr^ , (1) 

of the coupled Einstein and massless ghost scalar field equations. In ([1]) r denotes the ex- 
tra space-like coordinate, a > is a curvature scalar (we choice increasing warp factor), the 
determinant, 

V^=e^'^l'-|, (2) 

is M-independent and the oscillatory function u{t,r) has the form: 

u(t,r) = sin(a;t)Z(r) , 
Z(r) = Ce-^'^l'^lFa f-e-'^'^O • (3) 



Here u denotes the oscillation frequency, C is an integration constant and 1^2 (^) is the second- 
order Bessel function of the second kind. The solution ([T]) describes the brane located at the 
node of the standing bulk wave ([3]). In this paper we consider the case with one node, i.e. the 
function Y2{r) in ([3]) has a single zero at the position of the brane, r = 0. This is accomplished 
by the assumption: 

- ^ 3.38 . 4 

a 

In the equations of the matter fields, the oscillatory function ([3]) enters via exponentials. We 
suppose that the frequency u of standing waves is much larger than the frequencies associated 
with the energies of particles on the brane and in these equations we can perform the time 
averaging of oscillating exponents. Non- vanishing averages are expressed by the formula [T314T5] : 



(e^) = lo {Z) , (5) 

where Io{r) is the modified Bessel function. 

The brane at the node of the standing wave (the point where Z{r) and u(t, r) vanish) can 
be considered as a 4D 'island', where the matter particles are assumed to be bound. Indeed the 
system of 5D geodesic equations of motion for a classical particle (or a photon) reduces to [13J: 

r2 




+ a{l-V')\r\ = e , (6) 

where V is the velocity along the brane and the constant e corresponds to the energy of the 
particle per unit mass (for massive particles < e = E/m < 1 and for photons e = 1). From 
Q it is clear that the motion towards the extra dimension r is possible when 

e - a{l - V'^)\r\ > . (7) 

The term a(l — y^)|r| plays the role of the trapping gravitational potential and there exists 
maximal distance in the bulk |r|max ~ ^/a- In the standard brane approach with decreasing 
warp factor (a < 0) localization was achieved due to the fact that the extra space actually is 



finite [3lll]. In our case, the increasing of tlie brane warp factor, e^"''^', creates the potential well 
that confines particles. 

Let us now consider the localization problem of quantum fields with various spin using the 
brane solution ([1]) [T314T6] . To have localized fields on a brane 'coupling' constants appearing 
after integration of their Lagrangians over the extra coordinate, r, must be non-vanishing and 
finite. Using numerical solutions we shall demonstrate existence of normalizable zero modes of 
various fields on the brane within the model. 

2 Localization of bosons 

In this section, using numerical solutions, it is shown localization of massless modes of the 
spin-0, spin-1 and spin-2 particles on the brane. 

At first let us consider 5D massless scalar field. We assume that for the case uj ':^ E {E is 
the typical energy of scalar particles on the brane) variables can be separated [T3HT5], 



$ (x^) = 0(x^)^(r) 



Capital Latin letters numerate 5D coordinates, while Greek ones are used for 4D indices. Then 
the action integral for scalar fields can be split in two parts: 



S^ = -^Jd'x^ /^^^M^a^v^ = -iJd^^ \da(l)d''(l) I dre'^'^'^e - cf" j dre^^^^-^'^ 



, (9) 



where prime denotes the derivative with respect to extra coordinate r. We are looking for zero 
mode solution on the brane, ~ ^Wi^^" ^ qj^(\ assume that energy-momentum along the brane, 
Py, obeys the dispersion relation: 

e^-pI-pI-pI = (). (10) 

Than the 5D Klein-Gordon equation after time averaging reduces to [T314T5] : 

(e'"'\'k') -e^^\'\p\r)q = Q , (11) 

where 

P\r) = ((e-«) - l) {pI + pI) + ((e^«) - l) p^ (12) 

Close to the brane the function P'^{r), which describes oscillatory properties of standing 
waves, and the extra dimension part of wave function i;{r), behave as: 



P'ir] 






^Wl.^0 - e-"'^'- (13) 

So for the equation f ITT]) we use the boundary conditions: 

^(0) = 1 , ^'(0) = -4a . (14) 
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Figure 1: Scalar zero mode's profile in the bulk. 



Figure 1 shows the numerical solution of the equation ( ITTj) for these boundary conditions for 
the following values of the parameters: 

a = l, w = 3.38~a, p = O.OK a . (15) 

We see that i;{r), which governs of probability of scalar particles to reach the bulk, falls off 
from the brane to zero. 

The integrals over r in the action Q are convergent if integrand functions decrease stronger 
than 1/r. Figure 2 displays the products of r on these integrand functions. We see that they 
decrease, i.e. the integrals over r in the action ([9]) are convergent, hence 4D scalar fields are 
localized on the brane. 

Recall that the action of scalar fields (|9]) due to the determinant ([2]) contains increasing 
exponential factor (a > 0). This is the reason why in the original brane models [3ll5 the scalar 
field zero modes with the constant extra dimension parts can be localized on the brane only in 
the case of decreasing warp factor (i.e. a < 0). In our model the extra part of wave function 
( IT3|) is not constant, moreover, from Figure 1 it is clear that it strongly decreases and according 
to Figure 2 the integrals over r in the action ([9]) are convergent. 

It is known that the transverse traceless graviton modes obey the equations similar to the 
massless scalar fields in a curved background. Accordingly, the condition of localization of 
spin-2 graviton field is equivalent to that of scalar field considered above |13[[T^. 

In the case of 5D Abelian vector fields (the generalization to the case of non-Abelian gauge 
fields is straightforward) we are looking for the solution of field equations in the form [r ^ fT H ITB]: 

A(x^) = e-'''\''\g^pA^{xnv{r) , Ar{x^) = . (16) 

Than the 5D action of vector fields can be written as: 



Sa = -\J d'x^ g''''g''''J'MpJ'NR 



rWlrk^^H 



(r) \r\ g^°l^l 




Figure 2: Convergence of integrand functions in 5*$ 



-^ / d'x 



F^^F''^ I drv^ - 2A„A" / dre^'^^'^^v'^ 



(17) 



We require the existence of 4D vector zero mode on the brane, 



A^ {x") ~ e^e 



ipvX 



where p^ obey dispersion relation (jTOj) . Then 5D Maxwell equations after time averaging yields 
the single equation for v{r): 

(19) 
and we set the boundary 



[e^-\r\y'^ _p2^r)v = 



-'2a\r\ 



where P'^{r) is done by (IT^ . Close to the brane, v{r)\^,^Q r 
conditions: 

v{Q) = 1 , v'{0) = -2a . (20) 

From Figure 3 we see that the solution of the equation (IT9l) for these boundary conditions, or 
the probability of photon to leave the brane, falls off in the bulk down to zero. 

As for the case of scalar field the integrals over r in the vector functions action (fT7|) are 
convergent if integrand functions decrease stronger than 1/r. From the Figure 4 we see that 
the products of r with integrand functions decrease, i.e. the integral over r in (IT7|) is convergent 
and 4D photon is localized on the brane. 

In 5D brane models the implementation of pure gravitational trapping mechanism of vector 
field particles remains the most problematic |[9|I10]. The reason is that in the vector action (IT7|) 
the extra dimension parts of the determinant ([2]) and two metric tensors with upper indices 
cancel each other. Because of this in the original brane models the vector field zero modes 
(with the constant extra dimension part) cannot be localized on the brane for any sign of a. In 
our model the extra dimension part of vector field v{r) is rapidly decreasing function (Fig. 3) 
and the integral over r in (fT7|) is convergent. 




Figure 3: Vector zero mode's profile in the bulk. 



3 Localization of fermions 

Now we investigate the localization problem for massless fermions 
For Minkowskian 4x4 gamma matrices we use the Weyl basis. 



r 



/ 
/ 



7' 



-a' 
a* 



^^ ^ %^''^^^y^^ 



I 
-/ 



(21) 



and 5D gamma matrices T^ = /i^T^, where h\ = (g""!''!, e""'''!""/^, e""!*"!""/^, e^"l'''l+", 1). The 
inverse funfbein for the metric ([1]), we choice as: 



T'^ = /i^ t'^ , T" = i-f^ . 



(22) 



Non- vanishing components of the spin-connection in the background ([T]), which enter co- 
variant derivatives defined as: 



1 
_( 

4^ 



DA = dA + ^nfTsTc , 



(23) 



after time averaging are given by: 

(Q';) = a sgn{r) (h;) . (24) 

For the wavefunction of the bulk fermion field we use the chiral decomposition: 

^ (x^ r) = ^L {x") A(r) + ^r {x") p{r) , (25) 

where A(r) and p{r) are the extra dimension factors of the left and right fermion wavefunctions 
respectively. 
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vHr) \r\ 

i/'V) kk^"!'! 




Figure 4: Convergence of integrand functions in the vector field action ([T 



Then the 5D Dirac action for niassless fermions can be written as: 

+^RipL I dre^^^'^^p [A' + 2a sgn{r)\] - ^^t/jR f rfre^"l"lA [p' + 2a sgn{r)p] 



(26) 



We assume that 4D left and right Dirac spinors {'y^ipn/L = '^'4'r/l) are solutions of the free 
4D massless Dirac equations: 



R 



-ip^x" 








where p^ satisfy the dispersion relation (ITU|) and the constant 2-spinors L and R obey: 



-ip^x'^ 



{pt + aV*) L={pt- a'pi^ R = . 
Derived from (12^ 5D Dirac equation, after time averaging, takes the form: 

7*9t + (e"/2\ (^-d^ + ^^dy) + (e-") Yd,] ^ - 7^ [2a sgn{r) + 9,] * = 



te 



Using the expressions (l27j) and (|28|) it can be rewritten as the system: 

/ -e^l'^l [2a sgn{r) + dr] (r'Vi{r) \ ( p{r)R 

[ -a'V^ir) e''^''^ [2a sgn{r) + dr] ) [ A(r)L 



. 



(27) 



^28) 



(29) 



(30) 



Here we have introduced the functions VAr) 



n,(r) = ((e-"/2)-l)p, = [/o(Z/2)-l]p„ 
= ((e")-l)p, = [/o(Z)-l]p„ 



y 
VAr 



where Zir) is defined in (13]) 



A(r) 
P(r) 




(31) 



Figure 5: Profiles of the left and right fermion wavefunctions in the bulk. 



To find boundary conditions for the system (l30l) note that from the last two terms of (126|) 
we see that the action is even under the Z2 orbifold symmetry if A(r) is an even function of r 
and p{r) is odd (or vice versa) with respect to the transformation r -(r^ —r. We require, 



>^k{-r) = \k{r) , Pk{-r) = -pk{r) , 



(32) 



what can be accomplished if we impose that the bulk fermion wavefunction is even under 
five- dimensional parity [TTtfTS]: 



7^^(x'',-r) = -^(x",r) . 



(33) 



Then from the symmetry of the problem with respect to the transformation r ^ — r we get the 
following important boundary conditions: 



p(0) = , A(0) = 1 . 



(34) 



Solutions of the Dirac equations (l3Up with the boundary conditions (I34p are displayed on 
the Figure 5. We see that in our model the extra dimension part of the left spinor wave function 
A(r) has maximum on the brane and decreases in the bulk. While, as a consequence of ([3i 
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in our setup the right fermionic modes are absent on the brane and p{r) has maximum in the 
bulk outside the brane. 

Integrals over r in the spinor field action fl26p are convergent if integrand functions decrease 
stronger than 1/r. This feature for (p6|) is demonstrated on the Figure 6. So in our model left 
massless fermions are localized on the brane and right fermions are localized in the bulk. 



p(r) [A '(r)+2asign(r)' 
"AWle-*""-! 
■p^(r)|r|e^"i'i 
_A(r) [p'(r)+2a sign(r)' 

■p(r)lc^"'" 




Figure 6: Convergence of bulk integrals in Sij, ( l26l) . 



4 Massive modes 



In this section we want to estimate masses of Kaluza-Klein excitations. 

For simplicity let us consider the case of 5D scalar particles. For a massive mode the 
dispersion relation (ITU]) takes the form: 



77i2 2 2 2 2 



(35) 



where m denotes the mass of lightest KK mode on the brane. Then the equation for scalar 
particles (llip transforms to: 



/'e'^'^IH^'V _ g2a|r-| 



P^(r) -m' 



? = 



(36) 



It is more convenient to put 036p into the form of an analogue non-relativistic quantum me- 
chanical problem by making the change: 



For ilj{r) we find: 



^ r 



-2a\r\ 



ip{r) 



iIj" + e-'"l"l m' - V{r) ^ = , 



(37) 
(38) 



where the function, 



V{r) = 4ae-2"l^''l5(r) + Aa^e-^"^"^^ + P^ 



(39) 



is the analog of non-relativistic potential. Because this potential contains Dirac delta function, 
it becomes infinitely large at the point r = 0. This means that if we interpret ip{r) as a wave 
function, it must vanish at the point r = 0. Profile of the V{r) is shown on the Figure 7. 




■V(r) 



Figure 7: The effective bulk potential (l39l) . 

To show existence of mass gap between zero and massive modes we approximate V{r) by 
the rectangular potential and study quantum mechanical problem in potential well with the 
parameters: 

V(r) = { , ^ / ~ r( • (40) 

^ ^ I +CX) {r > L) ^ ' 

Here L is the width of the rectangular effective potential V{r), which can be estimated from 



Inside of the potential well the equation fl38|) obtains the form: 



< + e 



-2a|r| 



m^ - Aa^ 



V'i = 0. 



Thus in the first region, if 

m^ > 4a^ , 

the solution can be written in the form: 

^i(r) = Ci Jo (e-'^H/j) + C2Y0 (e-'^l^l/?) , (0 < r < L) 



(41) 
(42) 

(43) 



where Ci, C2 and Jq, Yq are integration constants and Bessel functions, respectively, and we 
have introduced the parameter: 

/3 = -Vm2-4a2 . (44) 
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In the second region, where the potential is infinitely large, we have: 

^2(r) = . (r > L) (45) 

The boundary and continuity conditions at the point of junction: 

z^i (0) = , V^i (L) = V^2 (i:) , i^[{L)=i^',{L), (46) 

as usual, lead to the quantization of the mass: 

MP) n(/3) 



Me-^(3) Yo{e-^(3) 
which can be satisfied if 



(47) 



/3 = -Vm2-4a2 = . (48) 

a 

So we find that in our model Kaluza-Klein excitations have discreet spectrum and the mass of 

the lightest KK massive mode is of order of the curvature scalar, 

m^2a . (49) 

This miens that KK states are very heavy since a correspond to the energy scale characterizing 
the brane as a topological defect in higher-dimensional space-time. 

5 Conclusions 

In this letter within the 5D standing wave braneworld model we numerically find profiles of 
massless matter field wave-functions in the bulk. Using numerical solution in action integrals 
we have demonstrated the existences of gravitational localization of all kinds of matter fields 
on the brane. For the case of fermions we found that, while left ones are localized on the brane, 
the right massless fermions are absent on the brane, they are localized in the bulk. We have 
also estimated masses of KK modes and found that, as in the standard single brane models, 
they are very heavy - of order of the energy scale associated with the brane. 
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